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We calculate spectrum of diffractively produced gluons in onium-heavy nucleus collisions 
at high energies. We show that it exhibits a characteristic dependence on nucleus atomic 
number A and energy/rapidity. We argue that this dependence offers a unique possibility 
for determining the low- a; structure of nuclear matter. Applications to RHIC, LHC and EIC 
experimental programs are discussed. 

I. INTRODUCTION 



Diffr active dissociation is one of the most interesting processes in high energy QCD. It played 

I2 j], je j, iy,y, 



a pivotal role in identifying early signatures of the gluon saturation 



16, 



17 



li], 3, 2q] 



10, 



21]]. While 



121 ] in deep inelastic scattering (DIS) at HERA Q, Q, Q, 
the measurements at HERA revealed the first indication that the gluon saturation has become an 
important effect, the possible measurements of diffractive dissociation in p(d)A collisions at RHIC 
and LHC as well as in DIS in the proposed EIC collider will be able to probe gluon densities deeply 
in the saturation region. This statement is supported by the recent phenomenological success of 
models based on gluon saturation that accurately describe the experimental data on total hadron 



multiplicities 



22 



23 



24 



36 



inclusive gluon production [26|, [27], 



37|, 



38, 



39 



40 



41 



28 



29, 



30 



31 



32 



33, 



34 



351 ] and 



431 ]. This motivated us to calculate multiplicity 



heavy quark production 

of diffractively produced gluons in coherent diffraction of onium on a heavy nucleus in a recent 



publication [44|. We observed that the diffractive gluon multiplicity is very sensitive to the low-x 
dynamics in onium. On the other hand, it showed only a weak dependence on the gluon density in 
the nucleus. The reason is that the total cross section is dominated by soft gluon momenta which 
are not sensitive to the short-distance structure on the nuclear color field. Consequently, in the 



present paper we set to calculate the diffracti ve g 
DIS has been discussed in many publications [14j, 



uon spectrum. Diffractive gluon 
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45|, 



46, 



47J, 



481, 



49 



50, 



51 



52, 



production in 



5J, 



5J. The 



2 



goal of our calculation is to calculate, for the first time, the diffractive gluon spectrum taking into 
account the low-x gluon evolution in all rapidity intervals, i.e. in the rapidity interval between the 
onium and the emitted gluon and between the emitted gluon and the nucleus. We believe that this 
calculation opens a new avenue towards the phenomenological applications in pA and eA collisions. 

The paper is organized as follows. In Sec. HJ we briefly review the result for the gluon spectrum 
given in our previous paper 441 ] . Eq. ([I]) represents the cross section in terms of the dipole density 



in onium n p (r, r', b, y) and the dipole-nucleus forward scattering amplitude N(r, b, y). In Sec. IIIII 
we deliberate about the behavior of n p (r,r',b, y) and iV(r,b, y) in various kinematic regions. We 
turn to analysis of the gluon spectrum for large dipoles r > 1/Q S in Sec. lIVI and small ones r < 1/Q S 
in Sec. [V] The results are summarized in Sec. IVII where we also discuss possible phenomenological 
applications. 



II. DIFFRACTIVE GLUON PRODUCTION IN ONIUM NUCLEUS COLLISIONS 




Consider a process of diffractive gluon production in onium-nucleus scattering such that the 
rapidity gap equals the produced gluon rapidity y. The corresponding fan diagram is displayed in 
Fig. [TJ In Ref. j^Jl we used the Mueller's dipole model 



of Kovchegov 



551 ] to generalize the quasi-classical result 



271 ] (derived independently in 53]) by including the quantum evolution effects. 
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The method is based on the principal idea of the dipole model that, due to the large difference 
between the coherence length of the low-x gluons in the onium light-cone "wave-function" and 
the nuclear size, we can split in the light-cone time the process of the low-x evolution in onium 
and the instantaneous interaction. Indeed, the coherence length of the low-x gluons is inversely 
proportional to x, whereas the size of the interaction region (in the nucleus rest frame) is 2Ra- 
In the large N c approximation the onium wave function decomposes into a system of independent 
color dipoles. Up to terms suppressed at low x, dipole transverse size does not change in a course 
of interaction with the nucleus. We therefore, are able to write the cross section for the gluon 
production as a convolution of the onium dipole density and dipole forward scattering amplitude. 
Let us introduce the following notations, see Fig. [2) transverse coordinates of quark, anti-quark, 
gluon in the amplitude and gluon in the c.c. amplitude are denoted by x, y, zi, Z2 respectively; 
gluon transverse momentum is denoted by k. In this notation the cross section takes the following 




FIG. 2: One of the diagrams contributing to the diffractive gluon production at the quasi-classical level. 
Notations are detailed in text. 



form 



da(k,y) a s Cp 1 
d 2 kdy = 7T 2 (2tt) ; 



J d 2 bd 2 B J d 2 r'ni(x-y,x! -y',B-b,Y -y) 



d 2 Zi e -ik-Z! 



Zi — X 



zi - y 



Zl — X' 



r\2 



m - y 



/\2 



[jV(x'-y',b,y) 



-iV(x' - zi,b, y) - N(y' -z 1 ,h,y) + N(psf - z u b, y)N(y' - Zl , b, y) 



, (1) 



where ni(r, r',B — b, Y — y) is the dipole density and N(r',h, y) is the forward dipole-nucleus 
scattering amplitude. Function ni(x — y, x' — y', B — b, Y — y) has the meaning of the number 
of dipoles of size x' — y' at rapidity Y — y and impact parameter b generated by evolution from 
the original 



equation 



56, 



dipole x — y having rapidity Y and impact parameter B 

57] 

<9ni(x-y,x'-y',b,y) _ a s N c f (x - y) 2 



5fJ. It satisfies the BFKL 



dy 2tt 2 J (x — z) 2 (y — z) 2 

[ni(x - z,x' - y',b,y) +m(y - z,x -y',b,y) - ni(x - y,x' - y',b,y)] , (2) 
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with the initial condition 

ni(r ) r',b,0) = tf(r-r')5(b) ) (3) 

where we denoted r = x — y and r' = x' — y'. 

he forward elastic dipole-nucleus scattering amplitude satisfies the nonlinear BK equation 



58, 



dN(x-y,b,y) a s N c f 2 (x - y) 2 

dy = ^J dZ (x-z) 2 (y-z) 2 ^ X " Z ' b < 

+N(y - z, b, y) - iV(x - y, b, y) - 7V(x - z, b, y)JV(y - z, b, y)] , (4) 



with the initial condition given by 



60] 



iV(r,b,0) = l- e -8 H %. (5) 
The gluon saturation scale is given by 

Q^ = ^Y^(b)^,l/r 2 ), (6) 

where p is the nuclear density, T(b) is the nuclear thickness function as a function of the impact 
parameter b. In the following we will assume for notational brevity that the nuclear profile is 
cylindrical. An explicit impact parameter dependence, which are required in the numerical analysis, 
can be easily restored in the final expressions. Accordingly, it is convenient to proceed by defining 
the quantity 

n p (r, r',y) = J d 2 b n p (r, r', b, y) (7) 

which satisfies the BFKL equation ([2]) with the initial condition 

n p (r,r',0) = 5(r - r') . (8) 

For the following calculations it is convenient to cast ([T|) in a different form extracting the 
explicit dependence on r'. First, we change the integration variable w = zi — y'. Then, introduce 
the following transverse vector 

X [iV(r', b, y) - JV(w - r', b, y) - N (w, b, y) + N (w - r', b, y)N(w, b, »)] . (9) 

Using Q, (HD can be rendered as 

da(k,y) _ a s C F 1 f , 2 , , 2 / , 2 , , uv , 2 



^7^2 | d 2 6d 2 B | dVm^^B-b^-^II^k^)! 2 , (10) 
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Now, contribution of the first term in the round brackets of ([9]) can be written as 

J d 2 w e-^ w ™~^ 2 [N(r', b, y) - N(w - r', b, y) - N(w, b, y) + N(w - r', b, y)N(w, b, y)] = 
- J d 2 w e *'(w-r') ?L [N(r', b, y) - N(w -r',h,y)- N(w, b, y) + N(w - r', b, y)N(w, b, y)] (11) 

where we changed the integration variable w — r' —* — w and used the fact that the amplitude 
depends only on the dipole size (and not on direction). Defining a new scalar function Q(r',k, y) 

as 

Q(r',k,y) = 

- J d 2 w e ikw -l [N(i>, b, y) - N(w - r', h, y) - N(w, b, y) + iV( w _ r ', b , y)N(w, b, y)] (12) 
and using (jlip we write © in the following form 

I(r', k, y) = -e- iW iV k Q(r', k, y) + iV k Q*(r>, k, y) . (13) 

Consequently, 

|I(r', k, y)\ 2 = 2|V k Q(r', k, y)\ 2 - e~^' (V k Q(r' , k, y)) 2 - e ik ' r '(V k Q V, k, y)f . (14) 

In the region where <5(r', k, y) is a real function, we can render (|14p as 

|I(r', k, y)| 2 = 4 sin 2 {^f^j (V k Q(r', k, y)) 2 , when Q(r', k, y) € (15) 

In terms of n p (r, r',y), (fT0|) reads 

dcr(fe,y) a s C F 1 /" 2 , . , 2 



5 A y dVn^r'^-y)!!^'^^)! 2 , (16) 



d 2 /tdy vr 2 (2vr) 2 

where Sa is the cross sectional area of the interaction region. At y = this expression reduces to 



the quasi-classical formula derived in 



48|, flEJ- 



So far we have been concentrating on a case in which the rapidity of the produced gluon y 
coincides with the rapidity gap Yq in a diffractive event. In this case the diffractive scattering 
amplitude iV/)(r, b, y, Yq) coincides with the square of the forward elastic scattering amplitude 
N(r, b, y). This case has the most phenomenological interest (since the invariant mass M of the 
produced system is dominated by "slow" gluons M 2 ~ k 2 /x). Still, at high enough luminosity a 
single hadron spectrum measurements should become possible. Therefore, a question may arise 
about the diffractive production of a gluon with y > Yq. Such process is shown in Fig. [3l In 
this case, the amplitude N(r, b, y) must be replaced by the off-forward diffractive dipole amplitude 



FIG. 3: Fan diagram for the diffractive gluon production in pA collisions with rapidity gap Y smaller than 
the gluon rapidity y. 



which explicitly depends on a coordinate of quark or anti-quark of parent dipole and the coordinates 
of the emitted gluon in the amplitude and in the c.c. one. 1 Investigation of properties of the off- 
forward diffractive amplitude would lead us astray of the main subject of this paper and hence will 
be discussed elsewhere. Let us only note here that in those cases when the coordinate of the gluon 
is the same on both sides of the cut, the off-forward diffractive amplitude reduces to the more 
familiar forward diffractive amplitude Nd(t, b, y; Yq). Since Nd(t, b, y; Yq) contains information 
about all possible pomeron cuts shown in Fig. [3] it can serve as a phenomenological model for the 
yet unknown off-forward diffractive amplitude. In this case the cross section for the diffractive 
gluon production takes form: 



dgP A {k,y) 
d 2 kdy 



oi s Cf 1 



7T 2 (2vr) 



^Sa / d 2 r' n p (r', Y-y) |I(r ; , k, y; Y )\ 2 



(17) 



where now in place of f|13f) and (I12j) we write 



I(r',k,y;y ) 




iV k Q(r', k, y- Y ) + iV k Q*(r', K V, *o) 



(18) 



1 In the case of inclusive gluon production, an off-forward amplitude was discussed in [6l| 
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and 



Q(r',k,y;Y ) 



dW k ' w - 



N*(r', b, y; lo) - ^ (w " r '> b > V5 ^o) 



-iV|(w, b, y; y ) + JV* (w - r', b, y; y )iV|(w, b, y; y ) 



(19) 



Amplitude N£>(r,h,y;Yo) equals the cross section of single diffractive dissociation of a dipole of 
transverse size r, rapidity y and impact parameter b on a target nucleus. It satisfies the Kovchegov- 



Levin evolution equation 



62] 



dN D (x-y,b,y;Y ) 



dy 

+ — ^— / a z 



(x - y) : 



27r<J(y-z)ln(|x-y|A) 



7T* 



(x - z) 2 (y - z) 2 

.2 

[N D (x - z, b, y; Y )N D (y - z, b, y; Y 



N D (pc- z,h, y;Y ) 



(x-yj 



(x - z) 2 (y - z) s 



4iV D (x - z, b, y; y )iV(y - z, b, y) + 2A(x - z, b, y)iV(y - z, b, y)] , 



(20) 



with the initial condition 



N D (r, b, y = y ; F ) = ^ 2 (r, b, Y ) . 



(21) 



Diffractive gluon production of the kind shown in Fig. [3] requires a dedicated study while in this 
paper we concentrate on the case y = Yq. 



III. DIPOLE EVOLUTION IN ONIUM AND NUCLEUS 
A. Dipole evolution in onium 

Dipole evolution in onium is encoded in the function n p (r, r',y) and is determined by solving 
the BFKL equation ([2]) with the initial condition ([3]). The result reads 

n p (r, r', y) = J dv (r/r') 1+2w , (22) 

where a s = a s N c /ir and 

x (y) = ^(1) - - iu) - \^{\ + w) , (23) 

with being the digamma function 

M-m- (24> 
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There are several cases when the integral (J22J) can be done analytically. Expansion near the 
maximum of x( u ) corresponds to the leading-logarithmic approximation. In this case we have 

n p (r,r',y) LLA « -±— /ZZ^ZZ e («J>-l)» e'^SM , a s y » ln 2 (r/r') , (25) 

where ctp — 1 = 4a s ln2. Alternatively, we can expand x{ v ) near one of its two symmetric poles 
at 2iu = dbl. This corresponds to the double logarithmic approximation depending on the relation 
between r and r'. The results for the dipole density read as follows 



n p (r,r',y) DLA 

and 

n p {r,r',y) DLA 



T l ( 2a ^) 1/4 2^/2a s yln(r'/r) I j / / / x ^ ^gx 

47r 3/2 r /4 ln 3/4 (r/ /r) e , r < r , ln(r /r) » a s y . (2b ) 



^J^^ e2 ^M^, r>/ , ln(r/r0 » asy . (27) 

B. Dipole evolution in a heavy nucleus 

In the region rQ s <C 1 where the forward elastic dipole-nucleus scattering amplitude N(r, b, Y) 
satisfies the BFKL equation it can be calculated similarly to the dipole density of the previous 
subsection. The initial condition in this case is specified by © expanded at small dipole sizes to 
the leading order. The result is 

N(v, b, y) LT = ^fj» e msXiu)V (rQso) 1+2i " - + ^ ■ ( 28 ) 

Analogously to the derivation of (|25h we obtain in the leading logarithmic approximation 

™- - ^T^k^ {¥) ^-*e-^ , a ,y » (-^-) , ( 29 ) 
and in the double logarithmic approximation 

ln 1 / 4 !^ 



r < 1/Q s0 , In — ^- > a s y . (30) 

Behavior of the scattering amplitude deeply in the saturation region rQ s 3> 1 can be found by 
noting, that with the logarithmic accuracy, the parent dipole r tends to split into two daughter 
dipoles w and r — w of different sizes: either u; ^ r ~ |w — r|or, symmetrically, |w — r| <C r ~ w. 



9 



Both give equal contribution to the integral over w. Restricting ourself to the case w <C r and 
doubling the integral we write the BK equation as follows: 



dN(r,b,y) a s C F 



dw 



[N(w,h,y)-N(w,b,y)N(r,b,y)]. 



(31) 



9y 7T Ji/qi w 

Now, for the reason that in the saturation region, the amplitude N(r, b, y) is close to unity we 



render (|3T|) as 

d{l-N(v,b,y)} a s C F 



dy 



TT 



l/Q 



2 w z 



{l-iV(r,b,y)} 



2o> s Cf , , 2 „ 2 



7T 



Hr 2 Qi){l-N(r,b,y)} 



(32) 



The saturation scale Q x {y) c an be found by equating the argument of the exponent in ([301) to a 



constant which yields 



63 



64] 



Q s {y) ~ Qsoe 



2a s y 



(33) 



63 



64 



65, 



Introducing a new scaling variable r = ln(r 2 Q^) we solve (]32|) and find the high energy limit of the 
forward scattering amplitude ([32]) 



N(r,b,y) = l-S e 



-r 2 /8 



661 ] (in the fixed coupling approximation). It reads 

1 



1 — Soe 



r > 



(34) 



where we approximated Cf ~ -/V c /2 in the large N c limit. So is the integration constant. It 
determines the value of the amplitude at the critical line r(y) = 1/Q s (y). 



IV. DIFFRACTIVE GLUON SPECTRUM: LARGE ONIUM r > 1/Q 8 

Now the stage is set for calculation of the diffractive gluon spectrum in various kinematic 
regions. Let us first analyze the differential gluon production cross section in the case of scattering 
of large onium r > 1/Q S . There are two interesting kinematic regions in this case depending on 
the relation between the gluon transverse momentum k and the saturation scale Q s . We consider 
these two cases separately. 



A. Hard gluons k > Q s 

To begin we need to calculate the function Q(r', k, y) given by (|12p . Note, that in the region 
w > 1/k the integrand is a rapidly fluctuation function. Therefore, the dominant contribution to 
Q(r',k, y) arises from dipole sizes w < 1/k. Consider now three possible cases, (i) r' < ^ < qj. 
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In this case splitting the integration region into two parts we write (|12p as 

r' d 2 w f l / k d 2 w 

Q(r',k,y)n —[l-N(r',b,y)}N(w,b,y)+ — 2JV(w,b,y). (35) 

JO w Jr 1 w 

In the second integral on the r.h.s. we neglected N(r', b, 0) as compared to iV(w, b, 0) since the 
amplitude is an increasing function of the dipole size and in most of the integration region w 3> r'. 
To determine the kinematic region that gives the largest contribution we note that when w <C <jj 
the amplitude scales as N(w,b,y) ~ w 2 . It follows, that the first integral in the r.h.s. of (|35p is 
of order r' 2 Q 2 whereas the second one is of order Q 2 /k 2 , i.e. the former is parametrically smaller 
than the latter. Thus, 

" 1/k d 2 w . _ . , 1 1 
— 2iV w,b,y , r' <-< — 
w z k Q s 



Q(r>, k, y) « / — 2JV(w, b, y) , /<-<—. (36) 



We obtain for the gradient 

V k Q(r',k,y) = -4n^N(k~ 1 k,b,y), r> < I < ±- (37) 

Eq. (|37p holds in the logarithmic approximation. Using (|15p we get 

|I(^k, y )| 2 =4^AT 2 (A ; - 1 k,b, y )sin 2 (^) , r'<i<i-. (38) 

In the second case (ii) -j- < r' < ^ and the third case (iii) | < Jr < r 1 there is only one 
significant integration region yielding 

' 1/fe d 2 w „, . , 1 , 1 ,1 1 



/••V" | 111 

Q(r',k,y)&[l-N(r',b,y)} jV(w,b,y), - < r' < — and - < — < r' . 

JO W K Lis K t/s 



(39) 



Therefore, 



V }i Q(r , ,k,y) = -27TjN(k- 1 k,b,y)[l-N(r',b,y)} , i < r' < -L and ~ < i- < / . (40) 

Substitution into (fT5|) yields 

|I(r', k, y)\ 2 = 4 JV^k, b, y) [1 - iV(r', b, y)] 2 sin 2 (^f \ , 

1 111 

- < r < — and - < — < r . (41) 

To calculate the differential cross section (|16j) we now need to integrate over all possible dipole 
sizes r' and orientations. Integration over the dipole orientations produces 
da a s C F c (4n) 2 u , / /"V* 



d 2 fc dy 



■«A^JV 3 (lb- 1 t,b > y) 



dr'r'n p (r, r', Y — y) (1 — Jo(& r')) 



1 /"VQ- 1 
+ 4 / fc ^Vn p (r,r', Y - y) (1 - Mkr')) j> . (42) 
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We restricted integration over r' to the region r' < 1/Q S since otherwise the integrand is strongly 

suppressed by [1 — JV(r', b, y)] 2 — ► 0, see (]4ip with (1341) or ([5j). To determine the largest contribution 

to the integral on the r.h.s. of (|42p we use the fact that the Bessel function Jq(x) ~ x -1 / 2 at x 3> 1 

and Jq(x) ~ 1 — x 2 /4 at x <C 1 and write the expression in the curly brackets as 

-i/fe i i rVQs 

dr'r'n p (r, r',Y — y) -k 2 r' 2 + - 
lo 4 4 J 1/k 



/■V* 1 1 /■!/<?. 

/ drV'n p (r,r',y-y)-A; 2 r ,2 + - / drVn p (r, r', F - y) . (43) 
Jo 4 4 J 1/k 



It follows from (|26p and ()27[) that n p (r, r', y) ~ r 2 / r ' 4 at r < r' and n p (r, r', y) ~ 1/r' 2 at r > r'. On 
that account, we determine that the first integral in (|43p is of order unity, whereas the second one is 
logarithmically enhanced by ln(k/Q s ) 1. A more accurate estimate is gained by substitution of 



(|27p and explicit integration over r'. Introducing a new integration variable C = 2^2(3^^ — y) In p- 
we have 

' 1/Q ' ,,,,,,, , v/2 & d( . v^f 



^ ' *V„„(r, r', y - „) „ ^ £ « e C = ^ [ erfi (V5) - «dl , (44) 

where erfi(z) is the imaginary error function defined as 

erfi(z) = — i eif (iz) , (45) 

Co = 2y / 2a s (Y — y) ln(rfc) and C* = 2^2^(1^ — y) ln(rQ s ). In compliance with the double log- 
arithmic approximation we must replace the imaginary error function by its asymptotic form at 
Co > 1 given by 

erfi(z) « -=- e z \ z » 1 . (46) 

Hence, keeping in mind that Co ^> C* we get 

rl/Qs i i , 

/ drVn p (r,r',r-y) « 2^ a (Y- y ) in(rfc) (47) 

A/k P 4tt3/ 2 (2a s (Y-y)ln(r£0) 1/4 V ' 

As expected, this integral is independent of Q s since the integrand is a steeply increasing function 

of -j. Finally, the cross section is procured by plugging ([4Tjl into (j4"2"|) 

^ r = ^^5 A iV 2 (A:- 1 k,b,y) * JV^-vJMrfc) r > J_ > 1 . 

d 2 fedy vr 5 /2 k 2 (2a s (Y - y) ln(rfc)) 1/4 Q s fc 

(48) 



B. Soft gluons fc < Q s 



We are now turning to analysis of soft gluon production by large onium. As in the case of 
hard gluons we wish to calculate Q(r',k.,y) in three different cases. First case corresponds to (i) 
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r' < jk- < jg, i.e. size of dipole emitting the triggered gluon is smaller than any other scale in the 
problem. We have 

r' d 2 w f 1/Qs d 2 w f 1/k d 2 w 

Q(r',k,y)«/ [l- N (v',h,y)]N( W ,b,y)+ / — 2 iV(w, b, y) + / — , (49) 

JO w Jr' w Jl/Qs W 

where we used the properties of the amplitude iV(w,b,y) as discussed after ([35]). The three 
integrals on the r.h.s. of (09]) is of order r' 2 Q 2 < 1, 1 and In & 3> 1 respectively. Evidently, the 
third one is dominating. Thus, 

Q(r',k,y) «2vrln^ (50) 

implying that 

In the second case (ii) J- < r' < ^ there are also three relevant regions of integration 

/•1/Q S d 2 f r> d 2 

Q(r',k,y) « / — [l-iV(r',b,y)]iV(w,b,y) + / — [1 - AT(r', b, y)] iV(w, b, y) 
JO w Jl/Qs w 

f 1/k d 2 w 

+ 1 —[l-N(r',h,y)]N( W ,b,y). (52) 



r' 



In the second and the third integral iV(w,b, y) ~ 1. The third integral is enhanced by In —n- and 
anyway it is the only integral that depends on k. Therefore, using f)34[) 

V k Q(r', k, y) « -2t3 [1 - iV(r', b, y)} = -2n* S e"S . (53) 

Consequently, 

|I(^,k,,)|' = l^Sje-^W) si „ 2 , _L < r , < 1 . (54 ) 

The third case corresponds to (iii) ^- < ^ < r'. There are now two relevant regions 

The second integral is enhanced by In jf- and, apart from the lower limit of integration, is the same 
as the third integral in (152p . Evidently, the k dependence of Q(r',k, y) is the same as in the case 
(ii), implying that (|54p holds in the case (iii) as well 



Q(r',k )2/ )« / ^[l-iV(r',b,y)]iV(w,b,y)+ / — [1 - N(r>, b, y)] N(w, b, y) . (55) 
Jo w Jl/Qs w 



W.KV)? = ^f- Sl e-l -*«*"> sin* , i- < I < r'. (56) 
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Essentially, what ([53]) and ([56]) tell us is that the region r' > ^ does not contribute to the 
cross section for diffractive production of soft gluon by large onium. Thus, the only contribution 
to the cross section stems from r' < ■A-. There are now two possibilities depending on the size r 
of the incident onium: (a) r > ^ > q- and (b) ^ > r > However, in both cases ^k- is the 
smallest size implying that the leading contribution to the cross section is the same in both cases. 
Expanding the argument of sinus in (15ip and substituting to (116p we have 

The dipole density is given by (I27p . Notice that since the largest contribution to the integral 
arises from dipoles of size r' ~ ^- (the integrand increase rapidly with r') we can approximate 
In jr ~ ln(rQ s ), neglecting contribution of very small dipole sizes r'. Thus 

da a s C F S A {2a s {Y - y)) l ' A 2y /2a B (Y- v )in(rQ s ) _ 1 . J_ /«s 

which holds for any relation between r and 1/fc. 

If we now wish to calculate the total cross section for diffractive gluon production at given 
rapidity y we have to integrate ()48|) and (i58j) over d 2 /c. Clearly, the leading contribution stems 
from the integral over soft gluons given by ()58|) . We attain 

(59) 



da _ a s C F (2a s {Y-y)) 1 ^ 2 j 2dts(Y ^ y)ln(rQs) lfo 
dy~8^ bA ln 3 / 4 (rg s ) ' VQ ' ' 



in complete agreement with the result obtained in our previous paper 



44]]. 



V. DIFFRACTIVE GLUON SPECTRUM: SMALL ONIUM r<l/Q s 

We now consider scattering of small onium on a heavy nucleus. We will again consider separately 
the two cases of hard and soft gluons. Calculation are facilitated a lot since we have already derived 
the function |I(r', k, y)\ 2 , which embodies information about the gluon emission and subsequent 
elastic scattering of the two intermediate dipoles w and r' — w off the nucleus. 

A. Hard gluons k > Q s 

Using (f38l) we obtain 

da _ a s C F 1 Q 4(4vr) 2 2 x 



d 2 kdy tt (27r) 5 



SA^-N^k^kMy) J dr'r'n p (r,r',Y-y)(l-J (kr')), (60) 
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where we integrated over orientation of the dipole r ; . To proceed we have to specify the relationship 
between the onium size r and the inverse gluon transverse momentum r- Assume that (a) r < 
k < -Q-. Then, integral over r' can be divided into the following four regions: (i) < r' < r, (ii) 
r < r' < ^, (iii) ^ < r' < -q- and (iv) J- < r' . To estimate the integral in each of this regions 
we use the same procedure as before (it is explained after (|43|) ). We find the following parametric 
dependence of the integral in these four regions: (i) k 2 r 2 , (ii) k 2 r 2 In j—, (iii) k 2 r 2 and (iv) r 2 Q 2 . 
Region (ii) gives the largest contribution. We have 

da _a s C F 1 c 4(4vr) 2 . .k 2 ' 



A (AirY ~ k l f 00 

S A ^LN 2 (k- 1 k,h,y) T dr'r' 3 n p (r,r',Y-y), (61) 



d 2 kdy 7T (27r) 

Upon substitution of (j26|) and changing to a new integration variable C, = 2\J 2a s (Y — y) In ^ we 
reduce the integral over r' to the imaginary error function as in (|44p . Following the same steps as 
those that led us to (14"T|) we derive 

f ' a ^S A r 2 N 2 (k-^,b,y) l - -eV 2 -^ 1 ^, r<i< ' 



d*kdy ^ * ' '*"(2a.(y-y)In ] J:) 1/4 ' * Q*' 

(62) 

Consider region (b) ^ < r < -q-. Repeating the same analysis as above we conclude that the 
dominant logarithmic contribution originates from the region x < r' < r. Accordingly, we use (j27H 
for the dipole density and neglect the Bessel function in (|60p . Doing the integral as explained in 
-(|47p we write 



J^T = ^SfS A l 2 N 2 (k^k,b,y) 1 eVWMrt), * < r < i_ . 

d 2 &dy vr 5 /2 A; 2 (2a s (y-y)ln(rA:)) 1/4 <2s 

(63) 



B. Soft gluons k < Q s 

In the case r < -q- < | formulas for |I(r',k, y)| 2 are given by (I5ip .(|54 p . (l56p . As was already 
mentioned, only small dipoles r' < 1/Q S contribute in the r 1 integral. We thus have two regions of 
integration: (i) < r' < r and (ii) r < r' < q-. The integral over the former is of order r 2 whereas 
over the latter it is of order r 2 In That being the case we derive 

da a s C F 1 [ l,Qs ,,, , , v . 4(2,) 2 k 2 r' 2 

Sa / dr r n p (r,r , Y - y) — -= — (64) 



d 2 kdy tv (2tt) 2 a J r py ' ' k 2 



^ (2^- y) ln^) 1/4 ' < «. <k {) 
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The total cross section is again dominated by soft gluons. Integrating (|64p over d 2 k such that 
k < Q s we find 

da o s Cf 



dy 47r 3 / 2 



Sa Oil r 2 



2a s (y-y)ln^ ; 



1/4 



l^/2a s (Y-y)ln^- s 



r < 



(66) 



again in agreement with our previous result 



44] 



VI. SUMMARY 



The differential cross section for diffractive gluon production is given by formulas (|48|) . (|58|) . 
(|62p . (|63p and (|65p . We can see that there are five distinct kinematic regions, which are really 
six. The behavior of gluon spectrum in these regions is sketched in Fig. [H To make the figure 

dN 1 



dN 1 dyfcdkJV„ 




FIG. 4: Sketch of diffractive gluon spectrum as a function of transverse momentum k in two cases r > Q~ 
and r < Qj 1 . Qualitative behavior in both regions is also indicated. xG(x — e v ~ Y ) is a "gluon distribution 
function" in onium, 7 is the anomalous dimension of the nuclear gluon distribution function and Nqo is a 
normalization constant. 



self-contained we indicated an approximate transverse momentum k and the saturation scale Q s 
dependence. 7 denotes the anomalous dimension of the nuclear gluon distribution. It varies from 
about unity at k 3> Q 2 /Q s o to 7 1/2 at Q s < k < Q 2 /Q s q- Fig. [J] teaches us that by varying 
the incident onium size with respect to the saturation scale we obtain different behavior of the 
gluon spectrum as a function of transverse momentum. In DIS the typical onium size can be 
varied by means of triggering on the events with different photon virtuality. Depending on the 
relation between k, Q s and r the gluon spectrum exhibits different pattern that allows a more 
direct measurement of the saturation scale Q s (y) (and hence the nuclear gluon density) than it 
is possible nowadays. The /c-dependence of hadron spectra is of course significantly modified by 
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the fragmentation process. On the other hand, dependence of hadron spectra on atomic number 
A is the same as for the gluon spectrum (it arises from the A-dependence of Q s , see ([6]), (l33j) ). 
The reason is that, as we explained in Introduction, the coherence length for the gluon production 
is much larger than the nucleus size, implying that the fragmentation process is independent of 
A. Consequently, A-dependence is a powerful tool in studying the nuclear gluon distribution. 



Likewise, energy/rapidity dependence is independent of details of fragmentation (see however [671 ]) 
and has been successfully used along with A-dependence for analysis of inclusive hadron production 
at RHIC. Therefore, energy /rapidity and atomic number dependence at different values of hadron 
transverse momenta allows access to information about the anomalous dimension 7, which is of 
crucial importance for understanding the transition region between the region of gluon saturation 
and the hard perturbative QCD. 

Similar arguments apply to the diffractive gluon production in pA collisions. In this case, 
however, there is a substantial uncertainty regarding the structure of the proton wave function. 
Diffractive gluon production in the case when the distance between the three pairs of valence quarks 
is about the same is strongly suppressed as compared to the case when the distance between one 
pair of quarks is much smaller than the distance between the other two pairs (quark - diquark 



configuration), see 



441 ] . In either case the A and energy dependence are given by Fig. U] (right or 



left panel). Since calculation of diffractive gluon production in pA collisions requires a substantial 
modeling of the proton wave function we intend to address it in a separate publication. 

To summarize, we calculated the spectrum of diffractively produced gluons at low- 2 in onium- 
heavy nucleus collisions. In the forthcoming publications we are going to apply our results for 
calculation of the diffractive gluon production in DIS and pA collisions. 
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